Delayed Scattering of Solitary Waves from Interfaces in a Granular Container 
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In granular media, the characterization of the behavior of solitary waves around interfaces is of 
importance in order to look for more applications of these systems. We study the behavior of solitary 
waves at both interfaces of a symmetric granular container, a class of systems that has received recent 
attention because it posses the feature of energy trapping. Hertzian contact is assumed. We have 
found that the scattering process is elastic at one interface, while at the other interface it is observed 
that the transmitted solitary wave has stopped its movement during a time that gets longer when 
the ratio between masses at the interfaces increases. The origin of this effect can be traced back to 
the phenomenon of gaps opening, recently observed experimentally. 
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The propagation of a perturbation in a chain of beads 
in Hertzian contact possesses soliton-like features, as 
first observed by Nesterenko 0. Several studies, ex- 
perimental 0, Q as well theoretical 0, @ have con- 
firmed the existence of such soliton-like pulses. De- 
spite the large amount of recent work on the subject 
S!!0I!HmEH[[lH0, the physics 
of granular media remains a challenge and new effects 
are there to be discovered and studied. Enlarging the 
number of (engineering) applications of such new effects 
needs a complete understanding of the dynamics of such 
granular media. 

The simplest granular systems are one-dimensional 
chains of elastic spheres. If they are in Hertzian contact, 
the spheres may be considered as point masses interact- 
ing through massless nonlinear springs with elastic force 
F = /ci5 3 / 2 , where 5 is the overlap of contacts and k is the 
spring constant (a function of the material properties) yij . 
Let Xi if) represents the displacement of the center of the 
i-th sphere from its initial equilibrium position, and as- 
sume that the i-th sphere, of mass m,, has neighbors of 
different radii (and/or mechanical properties). Then, in 
absence of load and in a frictionlcss medium, the equation 
of motion for the i-ih sphere reads 
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where it is understood that the brackets take the argu- 
ment value if they are positive and zero otherwise, ensur- 
ing that the spheres interact only when in contact. 

The interaction of a solitary wave with the boundary 
of two " sonic vacua" (meaning that the system does not 
support linear sound waves if not precompressed) was 
studied for the first time experimentally as well numeri- 
cally in [j| (see also 0,0,^1113 for a recent study). 

In this work we make a detailed numerical study of the 
propagation of solitary waves in a linear chain of beads 
composed of three "sonic vacua" , as shown in Fig. 1, that 
is, a granular container [2lll22| . This kind of systems are 
of interest because in them one can find the phenomenon 
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FIG. 1: Schematic granular container used in the calculations. 



of energy trapping. It will be assumed that all spheres 
have the same mechanical properties and that both ends 
of the chain are free to move. We have found that the 
scattering process is elastic at one interface, while at the 
other interface it is observed that the transmitted soli- 
tary waves take a long time to be released from one of 
the interfaces; this time gets longer when the ratio be- 
tween masses at the interfaces increases. The origin of 
this effect can be traced back to the phenomenon of gaps 
opening, recently observed experimentally 0E3. As 
far as we know, the effect found here has not been yet 
observed in experiments. 

Consider a set of spheres with two different radii a and 
b. It is known that adjacent spheres of radii a and b will 
interact with a force F = k a bS 3 ^ 2 , where 
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and E is the Young modulus and v the Poisson ratio of 
the bead material. 

We will consider the scattering of solitary waves in a 
system like the one of Fig. 1, consisting of a total of M 
beads. There are two set of beads with Ni beads located 
on the lhs, iV 2 on the rhs, both sets have beads with 
radii a and masses ra\. Between them there are L beads 
with radii b (a > b) and masses to 2 . Beads displacements 
are governed by a set of equations of motion that can be 
readily obtained from the successive application of Eq. 
||TJ, having in mind that the equation of motion for the 
first (resp. the last) sphere only includes the second (resp. 
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the first) term, in case when there is no wall (as we here 
assume). Spring constants are k^ in the middle, k aa at 
right and left hand sides and k a b at the interfaces of the 
granular system. 

In order to have realistic results, we shall assume that 
the system consists of stainless-steel beads (see Q for 
their properties), with radii a = 4 mm and 6 = 2 
mm. The number of beads is Ni = 30, N2 — 20 and 
L = 200. We also choose (3 = 10~ 5 m, 2.36 x 10~ 5 kg 
and a = 1.0102 x 10 -3 s as units of distance, mass and 
time, respectively. Through out the paper we assume 
that initially all beads are at rest, except for the first 
bead at the left side of the chain. This bead is supposed 
to have a nonzero value of velocity in order to gener- 
ate the soliton-like perturbation in the chain. We shall 
choose the following initial conditions 

Ui (0) = 0, i = l,...,M, in(0) = 101.02 0/a, 
71,(0) = 0, i = 2,...,M. 

This initial impact velocity corresponds to 1 m/s and 
therefore it is in the regime where plastic deformation 
can be neglected. The system is studied numerically 
by using an explicit Runge-Kutta method of 5th order 
based on the Dormand-Prince coefficients, with local ex- 
trapolation. As step size controller we have used the 
proportional-integral step control, which gives a smooth 
step size sequence. 



140 
120 
100 




50 100 150 200 250 
Bead number 

FIG. 2: Velocity of beads (in program units) as a function of 
bead number. Primary and secondary multipulse structures. 

As the solitary wave gets the interface a multipulse 
structure is generated but no backscattered solitary wave 
is observed. This last phenomenon has been explained 
by Nesterenko et al. [rjj as due to the opening of gaps 
in the vicinity of left interface. These effects originate 
from the discreteness of inertia and the nonlinearity of 
the interaction- they were first observed by Nesterenko 
and coworkers 

As the multipulse structure moves into the light sys- 
tem, there remains some energy behind the interface and 
after a while a second multipulse structure emerges, with 
similar characteristics than the first one but with less en- 
ergy [23|. This is shown in Fig. 2. (It can be shown that 



the opening of gaps in the vicinity of the left interface is 
also responsible for the emergence of this structure). 
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FIG. 3: Velocity of beads (in program units) as a function of 
bead number. Scattering of the leading pulse of the multipulse 
structure at the second interface, for different times: t = 1.303 
a (full line), t = 1.31 a (dashed line) and t = 1.32 a (long- 
dashed line). 

When the first multipulse structure interacts with the 
second interface it gets immediately scattered and both 
transmitted and reflected train of pulses appear in the 
right heavy system and in the light system, respectively. 
Part of this process is shown in Fig. 3. 

Now let us see what happens when the backscattered 
pulses interact with the left boundary of the granular 
potential well (see Fig. 4). It is observed that at time 
t = 2.32 a the leading backscattered pulse arrives at 
the interface. Contrary to what happens at the right in- 
terface, the transmitted pulse moves slightly forward till 
time t — 2.343 a where it "freezes". Even more intrigu- 
ing is the fact that only beyond time t = 2.462 a the 
transmitted pulse starts to move into the heavy system. 
In our original units of time this means that the transmit- 
ted pulse has stopped its movement during 1.2 x 10 -4 s 
approximately. This coincides with the fact that the sec- 
ond pulse of the multipulse structure approaches close 
to the interface. It is interesting to notice also that the 
backscattered pulses scatter without delay at this inter- 
face. 

If we allow for the first bead to have more energy at t = 
s, it is observed that the characteristics of the scattering 
process around the second interface is similar to the one 
observed before: the scattering is elastic. At the left 
interface the situation has not essentially changed; there 
still is a delayed transmitted pulse. 

To get more insight into this scattering process, let us 
observe the behavior of the velocity of beads as a function 
of time. We shall analyze that behavior at both interfaces 
to see the differences and try to find an explanation for 
this scattering. To that end we shall fix our attention on 
those beads around the interfaces. 

In Figure 5 we show the velocity of beads between 
beads 227 and 231 (i.e., around the right interface) in 
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FIG. 4: Velocity of beads (in program units) as a function of 
bead number. Scattering of the first pulse of the backscat- 
tered multipulse structure at the first interface, for different 
times: t = 2.335 a (dotted line), t = 2.35 a (full line), t = 2.36 
a (long-dashed line) and t = 2.462 a (short-dashed line). 
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FIG. 5: Velocity of beads (in program units) as a function 
of time (in units of a) for beads 227 (full line), 228 (short- 
dashed line), 229 (long-dashed line), 230 (dashed-dotted line) 
and 231 (dotted line). 



the time interval 1.303 a and 1.33 a. It can be straight- 
forwardly demonstrated that this behavior is analogous 
to that found in a system with only two "sonic vacua", 
in case a solitary wave, travelling in an unperturbed 
medium from the light to the heavy system, scatters from 
the interface. So it it important to stress that (in the case 
at hand) in the light system we have solitary waves trav- 
elling in it; they correspond (at least) to those leading 
pulses of the multipulse structure. 

In Figure 6 the velocity of beads 27 to 30 (i.e., around 
the left interface) is shown. It is notorious the big dif- 
ference between this behavior and the one found at the 
right interface. From here we deduce that the origin of 
the delayed behavior in the scattering process resides in 
the fact that, contrary to what happens at the right in- 
terface, around the left interface beads acquire a constant 
velocity in the interval t € (2.33a, 2.48a), where we have 
previously seen that the transmitted pulse has stopped 
its movement during a long time. It is also interesting to 
notice that bead 27 remains at rest during this interval. 

When changing the ratio of masses (in this case, the 
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FIG. 6: Velocity of beads (in program units) as a function of 
time (in units of a) for beads 27 (full line), 28 (short-dashed 
line), 29 (long-dashed line) and 30 (dotted line). 



ratio a : 6) one observes that the interval that takes the 
transmitted pulse to leave the interface increases. For 
example, one sees that by keeping the mass and initial 
velocity of the impacting bead as before, the intervals 
from the arrival to the left interface of the first backscat- 
tered pulse to the instant when the transmitted pulse 
starts to leave the interface are approximately 0.744 ms, 
for a = 6 mm and 6 — 2 mm, and 1.244 ms, for a = 8 
mm and 6 = 2 mm, respectively (compare with an inter- 
val of approximately 0.143 ms for the case a = 4 mm and 
6 = 2 mm. A numerical experiment with L = 50 beads in 
the interior of the container shows that this time interval 
does not depend on L). 
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FIG. 7: Velocity of beads (in program units) as a function of 
time (in units of a) for beads 27 (full line), 28 (short-dotted 
line), 29 (long-dashed line) and 30 (dashed-dotted line). 

In case a = 6 mm and 6 = 2 mm, Fig. 7 shows the 
behavior of the velocity of beads for beads 27 to 30 in the 
interval t £ (2.8, 3.392) a, where the transmitted pulse 
remains "frozen" at the interface. Notice that such be- 
havior for beads 27, 28 and 29 is essentially the same 
as the one observed in Fig. 6 for the same beads (in 
particular, bead 27 remains at rest). This confirms that 
the reason for the behavior observed in this work resides 
in the fact that some beads near the interface acquire a 
constant velocity during the interval where the transmit- 
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ted solitary wave retards its movement. Now, a constant 
velocity means that no forces are acting on them and 
therefore the phenomenon of gaps opening occurs. 

Using a detailed numerical approach, we have studied 
the scattering of solitary waves in a granular container 
consisting of three "sonic vacua" with Hertzian contact. 
We have found that the scattering process is elastic at 
the second interface, while at the first interface it is ob- 
served that the transmitted solitary wave has stopped its 
movement during a time that gets longer when the ra- 
tio between masses at the interfaces increases. At the 
same time, the reflected pulses appear to scatter elas- 
tically from the first interface. The opening of gaps in 
the vicinity of the left interface plays a crucial role in the 
observed behavior. The understanding of this kind of be- 
havior may be of help for applications of energy-trapping 
granular containers. 
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